Monovalent counterion distributions at highly charged water interfaces: 
Proton-transfer and Poisson-Boltzmann theory 
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Surface sensitive synchrotron-X-ray scattering studies reveal the distributions of monovalent ions 
next to highly charged interfaces. A lipid phosphate (dihexadecyl hydrogen-phosphate) was spread 
as a monolayer at the air-water interface, containing Csl at various concentrations. Using anoma- 
lous reflectivity off and at the I/3 Cs^ resonance, we provide, for the first time, spatial counterion 
distributions (Cs"*") next to the negatively charged interface over a wide range of ionic concentra- 
tions. We argue that at low salt concentrations and for pure water the enhanced concentration of 
hydroniums HsO^ at the interface leads to proton-transfer back to the phosphate group by a high 
contact-potential, whereas high salt concentrations lower the contact-potential resulting in proton- 
release and increased surface charge-density. The experimental ionic distributions are in excellent 
agreement with a renormalized-surface-charge Poisson-Boltzmann theory without fitting parameters 
or additional assumptions. 
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The electrostatics of aqueous solutions is a rich and 
fascinating topic full of unexpected and counterintu- 
itive phenomena that still presents noteworthy chal- 
lenges, both theoretically and experimentally, vital for 
a complete understanding of the physics of biologi- 
cal systems In recent years, there has been ample 
theoretical activity aimed at determining ion distribu- 
tions next to highly charged interfaces. It has been 
predicted that multivalent ions become strongly corre- 
lated next to the interface |^ Q, thus invalidating 
the traditional Poisson-Boltzmann (PB) theory Di- 
rect experimental verification for these correlations has 
been almost non-existent, and only recently is gradually 
emerging Although correlations among mono- 

valent ions are unimportant, the distribution of monova- 
lent ions at highly charged surfaces has also been contro- 
versial. In the context of hydration forces|0| recent theo- 
ries seem to favor interfacial restructuring of water pdlll^ 
leading to ion distributions that may significantly dif- 
fer from simple PB theory. First-principle predictions of 
surface-tension isotherms of surfactants assume the ex- 
istence of a relatively large Stern layer with a dielectric 
constant lower than that of pure waterll3l|. Furthermore, 
detailed theoretical and numerical studies have shown the 
importance of including other eflfectsQ- To settle these 
issues, it is essential to obtain precise experimental ion 
distributions, including distant points from the interface. 
This is also a necessary step for an unambiguous under- 
standing of those electrostatic effects that differentiate 
monovalent from multivalent ions. 

In this Letter, we report on experimentally deter- 
mined monovalent ion-distributions at highly charged in- 
terfaces, with an effective surface charge density in the 
CTo ~ 0.08 - 0.4^ range (molecular area of 40 - 180l^). 
Although our main interest in this study is in the con- 
text of biological physics, its relevance extends far be- 
yond that to basic aspects of intermolecular forces, elec- 
trochemistry and possibly to plasma physics. 



Poisson-Boltzmann (PB) yields ion distributions via 
two characteristic lengths, the Gouy-Chapmann length 
\gc ~ ksTe /inaQe, and the Debye screening length 

\d = \J s-Ke^ub ' "^^^sre T is temperature, ks is the Boltz- 
mann's constant, e is the static dielectric constant, and 
nf) is the bulk salt concentration. In the high charge 
limit, i.e., Xgc/^d << 1, PB theory predicts the dis- 
tributions next to the charged interface are practically 
independent of bulk concentration for at least the first 
10A(for 10-5 ^ nb ^ O.IM). 

Herein, we point to a well understood, but frequently 
overlooked issue concerning charged interfaces at aque- 
ous solutions. Nearly all biologically relevant molecules, 
including proteins, DNA, and many phospholipids be- 
come negatively charged by proton-release, the efficiency 
of which is given by a = 1/[1 + iqCp^^'-p^)]. Usually, 
the pKa < pH and almost all protons are dissociated for 
neutral pH ^ 7 {a ~ I). However, when such molecules 
form an interface, in particular a planar one, hydronium 
concentration becomes significantly higher than bulk at 
that interface, leading to a lower interfacial pH and to 
proton-transfer back to the interfacial molecules. Thus, 
the net surface-charge is reduced. Within PB the en- 
hancement is expressed quantitatively by the Boltzmann 
factor exp(— ^^), where "0(0) is the contact value po- 
tential. The effective charge at the interface is renormal- 
ized as follows HEI. 



1 + lOipK'^-pH)e t'BT 



(1) 



The potential at the interface, ■0(0), which can be influ- 
enced by ions in solution, is determined self consistently 
from the boundary-condition-equation sinh(0o/2) = 
— A£i/Aqp(0o), where 00 = eip(0)/kBT, and a renormal- 
ized Gouy-Chapman length = kBTe/2Trare, equiva- 
lent to the Grahame equation. The counterion distribu- 
tion is given by the Poisson-Boltzmann with a renormal- 
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ized Gouy-Chapman length (RPB). 

To experimentally determine the features of ion dis- 
tributions in water, we set up surface sensitive X-ray 
diffraction experiments from a well behaved and con- 
trolled Langmuir monolayer at the air/water interface. 
To extract the ion distributions, we employed the re- 
cently developed anomalous x-ray reflectivity technique 
for monolayers^. This basic type investigation has be- 
come feasible only with the advent of second generation 
synchrotron X-ray sources with novel insertion devices 
(i.e., undulator) and improved optics, which readily pro- 
duce variable-energy X-ray beams with brilliancies capa- 
ble of detecting a single atomic-layer even if not closely- 
packed. Another important advance in this regard is the 
development of liquid surface diffractometers first intro- 
duced by Als-Nielsen and PershanJ^. 

To manipulate ion bulk-concentrations, we used Csl 
(99.999%, Sigma Corp. Cat# 203033) solutions in ultra- 
pure water (experimental details handling monolayers for 
X-ray experiments are described in ll3|), taking advan- 
tage of the L3 resonance of Cs ions at 5.012 keV in 
anomalous reflectivity measurements. To control sur- 
face charge density, monolayers of dihexadecyl-hydrogen- 
phosphate (DHDP, see Fig. P) (C32H67O4P; MW = 
546.86, Sigma, Corp. Cat# D2631) were spread from 
3:1 chloroform/methanol solutions at the air-water inter- 
face in a thermostated Langmuir trough ITJ . DHDP was 
chosen for this study, since it forms a sim ple in-plane 
structure at high enough surface pressures [13 and its 
hydrogen-phosphate head-group (R-PO4H) has a pK^ 
= 2.1, presumably guaranteeing almost complete disso- 
ciation [P04]/[R-P04H] « 0.99999, with one electron- 
charge per molecule (ctq « 0.4 C/m^). 

X-ray reflectivity (XR) and grazing incident X-ray 
diffraction (GIXD) of monolayers at air/ water interfaces 
were conducted on the Ames Laboratory Liquid Surface 
Diffractometer at the Advanced Photon Source (beam- 
line 6ID-B, described elsewhere (l^J) to determine the 
structure of the monolayer and ion distribution^^ . 
Highly monochromatic X-ray beams (16.2 keV or 5.012 
keV; A = 0.765334 and 2.47374 A, respectively) were 
selected by a downstream Si double-crystal monochro- 
mator, and deflected onto the liquid surface to a de- 
sired angle of incidence with respect to the liquid sur- 
face by a second monochromator (Ge(220)and Ge(lll) 
crystals at 16.2 and 5.012 keV, respectively) located on 
the diffractometer. X-ray energy was calibrated with 
six different absorption edges to better than ±3 eV 
and subsequent energy scans at fixed Qz in the course 
of the present study, accurately confirmed the L3 en- 
ergy of Cs"*" [2l|. To extract the density profile across 
the interface from XR, a parameterized density profile 
p{z) = p'{z) + ip"{z) of the electron-density (ED) and 
the absorption-density (AD)(real and imaginary parts, 
respectively) is constructed by a sum of Error functions 

K-) = ^EErf(^)(p,-mi) + ^ (2) 



where N is the number of interfaces, pj = p'^ + ipj"] p'j 
and Pj" are the ED and AD of jth slab, Zj and aj are 
the position and roughness of jth interface respectively, 
Pn+1 is the ED of the solution. The AD profile is par- 
ticularly important at the Cs resonance (5.G12keV) as 
demonstrated below. The reflectivity is calculated using 
a recursive method [23| of the discretized density profile, 

Eq. ©.Him. 

Surface pressure versus molecular-area (tt — A) 
isotherms of DHDP at various Csl salt concentrations 
(rifc) shown in Fig.^ were used to control surface-charge 
density, in particular, to identify conditions under which 
DHDP surface-density (ctq = 1/^) is independent of rib. 
For TT > 0, the isotherm exhibits two distinct slopes, as- 
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FIG. 1: (color online) Surface pressure versus molecular 
area for DHDP for different Csl concentrations as indicated. 
Reflectivity and GIXD were performed at constant surface- 
pressures 40mN/m and 30mN/m. The dashed lines indicate 
the region in isotherm reported in the present manuscript. 

sociated with crytalline tilted and non-tilted acyl-chains 
with respect to the surface normal as identified by GIXD 
and rod-scans|2lj. In the present study we focus on 
the non-tilted crystalline phase (30 ^ tt < 40 mN/m), 
where the density variation at a fixed tt is less than 1.5%. 
This small variation in surface-density is corroborated 
by GIXD and rod-scan measurements of monolayers on 
Csl solutions in the 0.1-lO^^M range, which show two 
prominent inplane Bragg reflections consistent with the 
formation of 2D polycrystalline hexagonal svmmetrv|2l|. 

Figure 121 shows normahzed reflectivity curves, R/Rp 
(where Rp is the calculated reflectivity of an ideally flat 
water interface), for DHDP (tt — 40mN/m) on pure H2O, 
lO-^M, lO-^M, lO-^M, Csl measured at E=16.2keV. AU 
XR curves differ in the exact position and the sharp- 
ness of their minima, and the intensities of their max- 
ima. Similar reflectivity curves were obtained for tt = 30 
mN/m. The solid lines are the best-fit calculated reflec- 
tivities based on refined density profiles that show dif- 
ferences mainly at and below the phosphate head-group 
region. Since, the packing of DHDP is basically indepen- 
dent of salt concentration for tt = 40 mN/m, the reflec- 
tivity curves in Fig. |21 qualitatively show a strong depen- 
dence of ion distribution at the interface on bulk ion con- 
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FIG. 2; (color online) Measured normalizedXR curves and 
best fits (solid lines) for DHDP monolayers at various Csl con- 
centrations at surface pressure 40mN/m (curves are shifted, 
by a decade each, for clarity). 



centration, in agreement with RPB discussed above. To 
obtain counterion distributions, we combine the reflectiv- 
ities at and off resonance (5.012 and 16.2 keV) into one 
data set and refine structural parameters using a space- 
filling model and applying volume constraints of different 
constituents [3. l2ll l2a|. Figure 13 A) shows reflectivities 
of DHDP spread on IQ-^M at tt = 40 mN/m at 16.2 
and 5.012keV. The solid lines are calculated from gener- 
alized density p{z), obtained from parameter reflnement 
of a single model-structure for the combined data sets, as 
shown in Fig.OlB). The AD curve for 5.012 keV up to a 
normalization factor is practically the profile of the coun- 
terions at the interface (there is a minute contribution to 
the AD from phosphorous in the head group region). The 
difference between the ED's at and off resonance, normal- 
ized by Z[l — f'{Eres)] where, Z = 54 for Cs"*", gives the 
desired ionic distribution at the interface. Figure EIC) 
shows in a solid line the experimental Cs"*" distribution 
at the interface at 10~'^M. Similar distributions, at other 
bulk Csl concentrations, are shown (solid lines) in Fig. 
2JA). Using the space-filling model to analyze the X-ray 
reflectivities off resonance, and self consistently by in- 
tegrating the distributions (obtained by the anomalous 
reflectivity) along the Z-axis, the number of counterions 
per DHDP was determined. A compilation of the inte- 
grated number of ions at the interface is given in Fig. 
0^B) (square symbols), along with the calculated RPB 
(solid line). The values from PB theory are also shown 
in Fig. 2jB) (triangles connected with a dashed line). 

To account for the fact that PB equation assumes 
point-like charges and an ideally flat interface, we pro- 
pose to convolve the theoretical distribution n"^(z) with 
a Gaussian function as follows. 
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dz' 



(3) 



The convolution function is justified by the following as- 
sumptions: 1) interface fiuctuations are Gaussian and 2) 
the wavelength of the capillary waves at the interface 
are larger than molecular size. The width of the Gaus- 
sian F is determined by surface roughness, which is in- 
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FIG. 3: (color online) (A) Normalized X-ray reflectivities 
measured at 16.2 and 5.012 keV of DHDP monolayer spread 
on 10~^M Csl solution (7r=40mN/m). Solid lines are calcu- 
lated reflectivities using the ED and AD profiles shown in (B) . 
The two data sets were combined and reflned to a model with 
common structural adjustable parameters. (C) The solid line, 
obtained from the difference of the two ED's in (B), shows 
the experimental distribution of Cs"*" , the dashed and dashed- 
dotted lines are calculated by RPB and PB respectively con- 
voluted and non convoluted as indicated (non-convoluted cal- 
culations of PB and RPB are divided by 10). 



dependently determined from the reflectivity, and there- 
fore the convolution does not involve any new parame- 
ters. The calculated ionic distribution using RPB for 
— 10~'^M and its convolution are shown as dashed lines 
in Fig. injC) (F — 3.9 A is determined from our reflec- 
tivity experiments, and using pH-pKa = 2.5, within the 
range of uncertainty of the measured pH of our pure wa- 
ter), superimposed on our experimental data (solid line) 
with no adjustable parameters. Figure0lA) shows coun- 
terion distributions for three different n^'s (solid lines), 
with corresponding calculations of RPB and convoluted 
as in Eq. Q (dashed lines). The agreement between 
RPB theory and experiment with a simple smearing of 
the distribution, Eq. ijjjl is remarkably good, except for 
slight deviations at distances Z < —10 K away from the 
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FIG. 4: (color online) (A) Interfacial Cs^ distributions (solid 
lines and shaded areas) determined from anomalous reflectiv- 
ities as outlined in Fig. |21for various Csl bulk concentrations 
(shifted by 0.5M for clarity). Calculated distributions based 
on RPB are shown with dashed lines. (B) Number of Cs^ 
ions per lipid (« 4lA^) (square symbols). The dashed line 
is the PB integrated over the flrst IsA, and the solid line is 
obtained from RPB integrated over the same range. 



interface. PH. 

In this study, we have shown remarkable agreement 
of the RPB theory with experiment. We obtained both 
the distribution and the integrated number of monova- 
lent ions per charge at the interface over five orders of 
magnitude in ion-bulk concentrations. 

Our experimental results for the ion distribution are 
entirely consistent with water being described with a con- 
tinuum of bulk dielectric constant. Corrections due to 
finite size ionic radius, charge modulations, short-range 
interactions, image charges or water restructuring were 
not necessary for describing the experimental data, im- 
plying that such effects, if relevant, would change the 
distribution at distances shorter than ~ 3 A0|. The 
fact that the pK^ of the amphiphiles in our experiment 
is one of the lowest available also shows the dramatic ef- 
fects that the renormalization Eq. ^ (RPB) has at high 
surface charges. The understanding we gained with the 
monovalent ions is of critical importance to our ongoing 
investigations of charged interfaces at multivalcnt-ion so- 
lutions. 
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